

field = 'l028'

map = readfits('/d3/BGPS_V2/v2.0_ds2_'+field+'_13pca_map20.fits',hd)
lab = readfits('/d3/BGPS_V2/bolocat/release/v2.0_ds2_'+field+$
               '_13pca_labelmask.fits',lhd)




myps,'./grs_paper/bgps_example_field.eps',xsize=5

undefine,pr
pr=[-.12,0.75]
plotmap,map,hd,range=pr,ct=3,charsize=1.0,xr=[28.5,28.0],yr=[-0.2,0.3]
print,pr

plotmap,lab,lhd,OUTIMG=outlab,OUTHDR=outlhd,/noplot,$
        xr=[28.5,28.0],yr=[-0.2,0.3],XC=xc,YC=yc

labval = (outlab[uniq(outlab,sort(outlab))])[1:*]
nlab = n_elements(labval)
FOR ii=0,nlab-1 DO BEGIN
   
   cgContour,outlab EQ labval[ii],xc,yc,level=0.5,color='cyan',/over
   
   
ENDFOR
help,outlab







myps,/done



END
